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PREFACE 



This investigation was oonduoted at the United States Naval Post- 
graduate Sohool, Monterey, California as the thesis requirement for the 
degree of Master of Science in Aerology* 

For help and advice received in its preparation the author is indebted 
to Associate Professor G. J. Haltiner of the Postgraduate School. 
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Figure 1. 500 mb. Level Velocity Frofiiec Between 30° K and 

60° N along 120° W Meridian for 0300Z January 4, 1951 
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I. INTRODUCTION 



In the attempt to subject the state of the atmosphere to a mathe- 
matical analysis, and in particular, to explain the growth and behavior 
of wave like disturbances in the westerlies, the meteorologist is con- 
fronted with two major problems j 

(a) An imperfeot knowledge of the instantaneous state of the 
atmosphere and 

(b) The mathematical difficulties inherent in the attempt to 
integrate the equations of motion. 

Adequate observational data is available to permit us to draw some 
oonolusions as to the "average” conditions in the atmosphere, or, in the 
case of ohoosing a particular model to analyse, a "plausible" condition. 

The mathematical problem is simplified by considering the effect of 
a small perturbation on a known zonal flow and neglecting all differentials 
in the resulting equations of higher order than the first. This results in 
a set of linear differential equations — the perturbation equations. The 
mathematical details have, in mo3t oases, been further reduced by making 
oertain simplifying assumptions in the models. Although such assumptions 
introduoo artificialities when compared with the atmosphere, they essent- 
ially resolve the problem to the study of dynamic instability — that 
portion of the stability or instability contributed by the nature of the 
flow pattern itself. Kuo [4] states that, in studying two-dimensional 
nondiver gent flow. 
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Similar results would be obtained for the very general case 
with vertioal motion, solenoids, and divergence and convergence, 
if oertain approximations are made in reduoing the partial dif- 
ferential equations for a certain dependent variable to an ordi- 
nary differential equation. 

The technique of using the perturbation equations is not new, having 
been used over fifty years ago by Helmholtz, Rayleigh and others. In 1928 
Solberg showed that unstable waves can develop in a sloping surface of dis- 
continuity between two parallel streams of different density, supporting 
the results of his synoptic studies of wave development on the polar front. 

In 1939 Rossby developed his widely known equation for the speed of long 
waves in the westerlies in the special case of oonstant tonal motion in a 
homogeneous incompressible atmosphere. Sinoe that time, a great number of 
meteorologists have applied the theory to a variety of models. Maohta [sj , 
using Rossby’ s model, assumed a solution for the perturbation stream func- 
tion which allowed for a tilt to the trough line* Haurwitz also used Rossby’ s 
model but allowed for the curvature of the earth and obtained only slight 
modification to Rossby’s results. Charney £l] considered a special case of 
baroclinic three-dimensional flow* 

Kuo [4] developed general stability or iter ia for a two-dimensional, 
barotropic nondivergent atmosphere using general, symmetric velooity profiles* 
Haurwitz [3] considered pure shearing waves between two oonstant zonal 
two-dimensional wind fields on a rotating planet. The case of a shear zone 
of finite width was considered neglecting the variation of the ooriolis para- 
meter with latitude* This problem in considered briefly in this paper. 

It may be noted that in nearly all of the papers studied only a portion 
of the oomplete problem was investigated* Some are ooncerhed with the question 
of whether or not a perturbation will be unstable* Others have eliminated 
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potentially unstable waves from their solutions (as will be done in this 
paper) and proceed to establish the speed with which the long stable waves 
will move. This approach is entirely reasonable beoau3e if the wave is un- 
stable, the solution is valid only in the initial development and the sig- 
nificance of the real values of wave velooity is limited to this time. Also, 
as Charney £l] points out, the evaluation of the complex roots may present 
great difficulties. On the other hand, if the wave is essentially stable, 
the seoond order differentials whioh have been negleoted in the pertur- 
bation equations will oontinue to be small and the value of wave velocity 
obtains is valid for an appreciable length of time. 

In this paper, only the long, slow moving or retrograding waves are 
oonsidered. The frequency equation resulting from the solution of the 
boundary systems determined by a sinusoidal velocity profile and a linear 
shear tone are determined. The frequency equation for a probability curve 
profile was not fully developed due to mathematical difficulties encountered 
in the application of one of the boundary conditions. This matter is dis- 
cussed in some detail. 

The roots of the freauenoy equation in the case of the sinusoidal 
velooity profile are evaluated (for a particular wave length) for a synoptio 
situation in whioh the horitontal wind velooity profile was in fair agree- 
ment with the model. The results are found to be of the right order of 
magnitude and 6how that there may be more than one possible value of wave 
speed under a given 6et of conditions. 
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We shall now consider three velocity profiles and derive the equations 
relating wave speed and wave number. The perturbation equations to be used, 
as developed by Haurwitt [2] , assume incompressible flow. We shall further 

assume two-dimensional nondivergent flow. The variation of the coriolis 
parameter with latitude is assumed to be a linear funotion of y. (Appen- 
dix II). Under these conditions, the perturbation equations reduce to 

v/'(c~u) + v ( U" - 3 -H u c = O 

as derived in Appendix I* 

A complete analysis of this equation must include a study of the be- 
havior of v in the neighborhood of the singular point c - U. Such a point 
must occur if the wave is moving with a speed between U 0 and u max* Kuo 
has investigated this point for the case of general symmetric velocity profiles 
by including friotional terms in the above equation. These terms are small com- 
pared with the others exoept when o - U is small. 

He shows that if a singular point exists, a stable wave oan exist only 
if U"-3 = O at some point on the velocity profile. If this condition is 
met, the speed of such a stable wave is the value of U at that point. Kuo 
calls this the critical velocity. Any wave moving with less than this 
velocity (but greater than U 0 ) will be unstable, while those moving with a 
velooity greater than the oritioal velocity will be damped. He further shows 
that for the long slow moving waves (c <U Q ), the amplifying or dampening 
factors are small or tero and such waves are essentially neutral. 
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In the solutions to be disoussed, no olosod form was found and the 
power series solutions obtained are valid only for (o - U) 0, the slow 
moving or retrograding waves.* 

* Kuo has shown that the minimum wave length of such waves is about 

5200 km under typioal shear oondltions* 
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II. SINUSOIDAL VELOCITY PROFILE 



The first problem to be considered is the solution of equation (l) 



where 

U = U. +0, a xu Xr ^ o 4 ^ 

0 - 0 o 

^ Z lr 

A series solution for v in the sinusoidal region is obtained in the con- 
ventional manner as shown in Appendix III. It involves two arbitrary 
constants whioh must be eliminated by the imposition of two boundary 
conditions * 

Euo [4} assumes that for symmetric velooity profiles and symmetric 
boundary conditions, there will be either symmetry or antisymmetry in the 
amplitude of the perturbation stream function, and therefore only half the 
*one need be considered. He then points out that it would be unlikely to 
have v * 0 along the line of maximum wind, thus reduoing the study to one 
of symmetric disturbances. 

Kuo's development is in terms of , the perturbation stream funotion 
instead of v. It oan be shown that, under the assumptions involved, equa- 
tion (l) in either variable is identical. 

The above statements result in the conclusion that v has a maximum 
where U is a maximum; that is. 




4 6 



( 2 ) 



(In this and subsequent work, the quantity in parentheses refers to the 
point at which the funotion is to be evaluated) 
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This equality is used as the first boundary condition and results in 
the elimination of one of the arbitrary constants. 

Where y < 0, in the field of constant tonal wind, it can be shown 
that the solution of (1) is 

V = K e 

Then , 

V - n. V = o 

This condition must be true at all points in the tone, therefore it is 
true at y ~ 0. Since y * 0 is common to both this tone and the sinusoidal 
tone, the statement is true at this point when applied to the sinusoidal 
2 one. The seoond boundary condition therefore is 

V ‘ ( o') - A. V (o) = o 

The imposition of the boundary conditions is shown in Appendix III 
and results in a fourth degree equation in d, where d s c - U 0 « The four 
roots of this equation may be evaluated, by rather tedious methods, after 
selection of the physical parameters involved. 

Two synoptic situations in which the velocity profiles at the 600 mb. 
level were in fairly good agreement with that used in this development are 
considered as examples. These are illustrated on the following page. The 
parameters of the upper profile were evaluated and the equation determined 
for a wave length of 5000 km. This resulted in 

4 s - H S' JL - 1 8 1 0 4 * - I* 2 o 4. - 5loo - o 
The roots of this equation are (approximately) 

d * -2 m/soo., -1.8 m/seo., -22 m/sec., 71.5 ny/seo. 
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The last root may be discarded sinoe we are not considering waves whioh 



move faster than minimum wind velocity. This gives 
of wave speed as 13.4 m/sec., 13.6 m/see., and -6.6 
stable since all the roots are real. 



the possible values 
ry^seo. The wave is 
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III. PROBABILITY CURVE VELOCITY PROFILE 



The second velocity field to be considered in the solution of (1) 



is 



-Ay" 

u = +d, e 



( 3 ) 



By varying the parameters, this profile may be adjusted to give a good fit 
to a large number of symmetric profiles in the vicinity of the jet. Tho 
differential equation is solved by the series method in Appendix IV. The 
first boundary condition is the same as in the previous problem and when 
applied results in a solution for v involving one arbitrary constant. 

The seoond boundary condition, as before, is 



v* - rv s 0. 



However, the ’’point" at whioh this must be applied is at infinity. The 

use of this condition therefore requires a transformation of the equation 

in which y : i , The point at infinity is then transformed to the origin 
s 

and the boundary condition may bo applied. However, the mathematics here 
is rather complicated and the actual expression for c in terms of c< was 
not obtained. 
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IV. LINEAR SHEAR ZONE VELOCITY PROFILE 



The final problem is oonoerned with the solution of (1) for a wind 
profile which is defined by the equations 



This model differs from the previous two and from those considered by Kuo 
in that it is not symmetrical. It is also obviously less realistic. The 
solution for v in the shear zone is given in Appendix V. 

By similar reasoning as in the proceeding problems, the boundary con- 
ditions are found to be 

v'fo) - /l,V(o) =. o 



where r-^ refers to the zone of velooity U^j and r-^ to the zone of Uji* 

By applying these boundary conditions as indicated in Appendix V, a 
rather complicated equation containing an implicit relationship between o 
and results. The difficulties of finding the roots of thi6 type equation 
are so great that no attempt was made to obtain numerical results for a par- 
ticular model. 
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V. CONCLUSIONS 



As with any second order differential equation with variable coef- 
ficients, the equation solved in this paper does not, in general, yield 
a solution in terms of elementary functions. It is quite possible that 
the problems considered here do possess such solutions but none was found. 
An attempt was made in each oase to transform the equation so as to obtain 
a solution in terms of Bessel or Legendre functions, with no sucoess. An 
arbitrary selection of wind profile will, in general, be met by the same 
difficulties. 

The perturbation equations are, therefore, readily integrable only 
if the velocity profile is chosen specifically to achieve this purpose. 

A particular but arbitrary profile such as one of those ohosen for this 
paper will, in general, yield an equation in which the algebraio diffi- 
culties in obtaining numerical results are of such magnitude as to make 
them praotically worthless. These difficulties would be increased in the 
case of complex roots. 

The restrictive assumptions made in all published studies in this 
field have been justified, to a certain extent, but must not be forgotten 
when attempting to apply results to the atmosphere. The prooessess of the 
real atmosphere are much more complex than we have assumed. The three 
dimensional nature of the atmosphere, divergence, vertical motion, and 
non -homogeneity certainly play an important role, and often a dominant 
one in the development and propagation of disturbances in the zonal flow. 
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APPENDIX I 



DERIVATION OF SECOND ORDER ORDINARY DIFFERENTIAL EQUATION 
IN v AND y FRCM PERTURBATION EQUATIONS 



The perturbation equations for two-dimensional nondivergent flow 
in an incompressible fluid where the undisturbed flow is zonal are: 



3 k . , . ^ k . 

IT * u 71 
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+• M 


-»«■ -V 7* 

\ <3 x 


4* * U 4? 

dt d X 
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4- Jr H. - - p — 

f 3-J 


d U . <} v 






dx a 3 


— u 





(5) 

( 6 ) 

(7) 



In order to study the velocity of wave perturbations in the x direction. 



we may assume the following form of solutions 

c « (y-c * ) 

- u.( tj) Si 
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L <( r-e-fc) 

V s - v c y) a 



(r-cr) 
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( 10 ) 



Then from (7) and (8), 
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Substituting (11) and (10) in (5) we obtain 
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Talcing the partial derivative of (12) with respect to y find equating it 

to the left hand member of (6), using (9) and (11), we get an equality in 

L o< cx-cr) 

which each term contains the coefficient J) . Dividing through 

by this quantity reduces the equation to one in which the only dependent 
variable is the amplitude factor of v, v(y). For sake of brevity, v(y) 
will now be designated as simply v. The differential equation then 
reduces to 

V "(C- U) + = 6 

This may be considered as an ordinary differential equation, but it must 
be remembered that the final solution of v contains the periodicity factor 

14 (X -C tr} 
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APPENDIX II 



VARIATION OF CORIOLIS PARAMETER WITH LATITUDE 



Let the ooriolis parameter he represented by 

5 = 5 .+ *} 

where J is assumed constant in the zone under consideration, 

hut 5 J Z O 

and - Z <ju co-6 <d 

From the sketch we see that 




Substituting this value of ^ in (A) we have 
— ; J- Z CO <D 

ft 



Then B 
(A) 

0 
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APPENDIX III 



DEVELOPMENT OF FREQUENCY EQUATION FOR SINUSOIDAL VELOCITY PROFILE 



We shall now solve the equation obtained in Appendix for the case 



of a sinusoidal velooity profile: 
L) - ID 0 4 (_) , l r 'j 



0 * J * V 



(13) 



The sinusoidal profile is bounded on both sides by a field of constant 
velooity, U Q . 

Substituting (13) in (l), we have the following equation to be 
solved : 

u" ( cl - u, i.y) ^ ll) - 0 (14) 

where 



d - C - 0 
<C - 

b _ B - < " A. 

Expressing sin by in series form 

ju. i./ .. iy _ W . IhL _ 



7 -■ *7 



i I 



a'i 



(15) 



and assuming a power series solution for v, 

>/=&*> +- a , y ^ y t- - - ~ ' 



( 16 ) 



we may substitute (15) and (IS) in ( 14 ) and obtain an equation in positive 
powers of y equal to zero. Upon equating the coefficients of each power 
of y to zero, we find that 
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This expression for v involves two arbitrary constants whioh may be elim- 
inated by the imposition of two boundary conditions* The first condition 
to be applied, as shown on page 6, is that v have a maximum where U is a 
a maximum, or. 
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- 0 



Taking the first derivative of (17) with respect to y, setting 
and equating the result to zero, we have 
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(18) 



and upon solving for aj we see that 

q <t h t- J + 0 t h- tt 1 * ( h +- o J ) 

Q> 1 InTF - w -n-" 



(19) 



The resulting solution for v, neglecting, for the moment, the third and 
higher powers of y, is 



s/ -- Ol. 
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(20) 



3y imposing the aeoond boundary condition (see page 7), 
v ( o) - -' 11 /( 0 ) — 0 

we obtain 






MTT^-t -4 t-U, ^ 

8 -t'~ A" - V. (i TT x 




( 21 ) 



This step shows that only the first two terms of the series solution for v 
are significant in the development of the frequency equation so the powers 
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of y higher than the first could have been ignored after having deter- 
mined equation (19). However, it should be noted that the solution of 
a^ in equation (1?) is approximate since only the first four terms of 
the general solution for v were used in its determination# It follows 
that the frequency equation is approximate# 

Squaring equation (21) and collecting the coefficients of powers 
of d we have 
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APPENDIX IV 



DEVELOPMENT OF FREQUENCY EQUATION FOR PROBABILITY 
CURVE VELOCITY PROFILE 



In this section, the differential equation 

v ( C - u ) + C 04 u ^ " 0 (1) 

is to be solved for the case in •which the undisturbed wind velocity 
profile is defined by the probability ourve 

0 = + U, A (23) 



A power series eolution for v is assumed. 



v - a „ + 4, y < a t y 



-*/ 



(24) 



and £ is expressed in series form, 

-v\ v kY 

& - / - ft / + 3 / 



(25) 



The series expressions for v and U are then substituted in (1). The 

coefficients of each power of y are equated to zero and we may now obtain 

expressions for any a n in terms of a Q and a^s 

- ZkU, - 0 

^ ^ O ~~ 

2- - O 



- 4. 



- au,-^c. 6 

l> ( c ; 



where 0,-t- U, and 



a 0 and a^ are arbitrary constants. The general solution for v is therefore 

(26) 



, _ „ , , « ^m-Uky, -C< l 'C-0 «- *-/. ,,, 

yj - a u 1 -u, y +. 4„ ! a ikb, -w'-c.- 3 } 

i(o^- C ) 1 +*■, -77^ : / + 



ls>(0^ . c ) 
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Vie may eliminate one of the constants by the imposition of the boundary 
condition 

\/ ' ( 0 ) - O 

The use of this condition for symmetrical velocity profiles is discussed 
on page 6# 

Taking the first derivative of (26) with respect to y, and setting y = 0, 
we obtain 

CL, = o 



Equation (2o) then reduces to 

- zk u. -Cq< L - 9_ 






s *4 



i 4 



2. ( Urn -C) 



y 



4- 



As mentioned on page 10, the application of the second boundary 
condition, 

~ vu v - 0 



leads to mathematical difficulties and the elimination of a^ and consequent 
development of the frequency equation was not accomplished. 
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APPENDIX v 



DEVELOPMENT OF FREQUENCY EQUATION FOR A LINEAR SHEAR ZCNE 



In this seotion the differential equation 

( c - o ) + \/ ( </' - 3 - 



( 1 ) 



is solved for the case in which the undisturbed wind velooity profile is 
a linear shear zone. The arbitrary constants are eliminated and the fre- 
quency equation developed for the oase in whioh the shear tone is termi- 
nated, on eaoh side by a field of constant zonal wind. 

The shear zone is defined by 
l> — . + )< / i^heK k - — — — 

*■ (27) 



A power series solution for v is assumed 
V- a o + a , y -i- a t y - 

and this series is substituted in (1). The coefficients of eaoh power 

of y are equated to zero and the following relationships obtained* 

<* 1 + b 



A, G* 



kb d U L A ^ B 

* A. , - 



’ Ud 

The series solution for v is then 



where J. 1 C - 0 ( 



a 0 + a., y +. a e 



■Ji- 0 



/ 



UJ 1 - 



d + 1 J 



W 1 



i +■ • • • (28) 



As explained on page 11, the boundary conditions to be imposed are 
v(o)-A.,V(o) = o 



\/ 



(£) f A., V(j) - 0 
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By substituting (28) in (29) and letting y - 0 we find that 



CL , — i A o 

Equation (28) then reduoes to 



V = a« 



+ A. / 



c< B 
Z <4 



1 



K 3 + A,d v/ 3 

‘ 



b d L 



( 31 ) 



To impose the second boundary condition, we substitute (31) in (30) and 
set r* . After dividing through by the common faotor a 0 , we obtain the 
frequency equation whioh beoomes, upon rearranging: 

Is, tl'fl,* A..+ A./v,, Jf) -h (, 4? U^-4 fi) 

3 *„ J ’'(e^ cJ + 5^ + 3 J 1 '[k 5 -t h, 

V P ** [ K& + A-, C* ^ — ° ' 



This is approximate since fourth and higher powers of y have been neglected. 
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